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ABSTRACT 

Using the Hubble Ultra Deep Field Near Infrared Camera and Multi-Object Spectrometer 
(HUDF-NICMOS) UDF05 parallel fields, we cross-matched 301 out of 630 galaxies with the ACS 
filters V 0 O 6 and 2850 , NICMOS filters Jno and H 160 , and Spitzer IRAC filters at 3.6, 4.5, 5.8 , and 
8.0 fan. We modeled the spectral energy distributions (SEDs) to estimate: photometric redshifts, 
dust extinction, stellar mass, bolometric luminosity, starburst age and metallicity. To validate 
the photometric redshifts, comparisons with 16 spectroscopic redshifts give 75% within A < 0.2, 
which agrees with the sensitivities expected from the Balmer-break in our dataset. Five parallel 
fields observed by NICMOS have sensitivities in the Higo-band of 80% at m,\B = 25.4 and 50% 
at rriAB = 26.7. Because the sample is Hi6o-band selected, it is sensitive to stellar mass rather 
than UV luminosities. We also use Monte Carlo simulations to determine that the parameters 
from the best-fit SEDs are robust for the redshift ranges 2 > 1.3. Based on the robustness of the 
photometric redshifts, we analyze a subsample of the 301 galaxies at 1.3 < 2 < 2 (35 objects) and 
3 < 2 < 4 (31 objects) and determine that Lb 0 i and the star formation rate increase significantly 
from 2 ~ 1.5 to 4. The Balmer decrement is indicative of more evolved galaxies, and at high 
redshifts, they serve as records of some of the first galaxies. Therefore, the galaxies in this sample 
are great candidates for future surveys with the James Webb Space Telescope and Atacama Large 
Millimeter Array. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: 
photometry 


1. Introduction 

As deeper infrared surveys have become avail- 
able, the ability to detect fainter galaxies with 
prominent spectral features in the rest-frame op- 
tical is proving an important key in the evolution 
of the early universe to the present epoch. Mul- 
tiple wavelength photometry (optical-infrared) is 
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often used to determine the properties of high red- 
shift galaxies, such as star formation rates, stellar 
mass and bolometric luminosity from inferences 
of nearby observations. The physical properties 
of nearby galaxies fit along the Hubble Sequence 
when observed in optical wavelengths. These 
properties often correlate with gas and dust con- 
tent, star formation rates, ages and luminosities. 
At around 2 ~ 1 this classification scheme quickly 
fails. The relevance of using the known properties 
of nearby galaxies to compare with high redshift 
galaxies has been revisited many times showing 
that certain characteristics may be compared di- 
rectly (e.g., Weedman & Huenemoerder 1985; Hi- 
bbard k Vaeca 1997; Heckman et al. 2005; Lotz 
et al. 2006; Overzier et al. 2008; Petty et al. 
2009; Overzier et al. 2010). 
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One method of combining the light at multi- 
ple wavelengths to determine the integrated stellar 
properties of galaxies includes the fitting of photo- 
metric data with the overall shape of galaxy spec- 
tra, or spectral energy distribution (SED). The 
technique is well-established and first proposed by 
Baum (1962) for obtaining photometric redshifts 
of cluster galaxies. Since the mid-90’s photometry 
became trusted as a reliable method for determin- 
ing redshifts (e.g. Mobasher et al. 1996; Lanzetta, 
Fernandez-Soto, & Yahil 1998; Hogget al. 1998). 
It has proven to be invaluable, especially for deep 
and large surveys such as the HDF North and 
South, since spectroscopy is time consuming and 
not feasible for large galaxy surveys. Photomet- 
ric redshifts are derived from strong spectral fea- 
tures, such as the Lyman-break (912 A), the 4000 
A break, or the Balmer-break (3646 A). 

One main difficulty with observations of high 
redshift galaxies is correlating the observed prop- 
erties with the appropriate rest-frame. Infrared 
(IR) observations of galaxies show the rest-frame 
near-IR (NIR) and optical for z > 1, while optical 
observations show the rest-frame ultraviolet (UV). 
At high redshifts, observing the rest-frame IR to 
identify thermal emission from star-forming re- 
gions becomes problematic in the FIR, because of 
the lower resolution images. Although optical ob- 
servations provide better resolution images, dust 
significantly masks the light at rest-frame wave- 
lengths blue-ward of ~ 3000 A . This significantly 
limits the types of galaxies resolved at high red- 
shift where photometry may be analyzed. 

SEDs are determined from templates of local 
galaxies and this presents a problem when at- 
tempting to fit models for higher redshift objects. 
Despite these complications, many important dis- 
coveries have been made concerning the nature of 
high redshift galaxies using SEDs. One of the most 
observed are Lyman-break galaxies (LBGs), which 
are compact, starbursts that are rest-frame UV se- 
lected (see Steidel et al. 1996) for the presence 
of a prominent Lyman-break at 912 A. The break 
is created by O and B stars that ionize hydrogen 
from the ground state. They have been identified 
from s ~ ltoz > 5 (e.g., Steidel et al. 1996; 
Bouwens et al. 2006; Burgarella et al. 2007). 
Older stellar population galaxies, such as Balmer- 
break galaxies (BBGs) have a drop in flux at the 
Balmer 3646 A. line. The Balmer-break marks the 


termination of the hydrogen Balmer series and is 
strongest in A-type stars, making it sensitive to 
stellar age. 

In the optical rest-frame, the Balmer-break is 
often used to select for galaxies with more evolved 
stellar populations (e.g. Baum 1962; Franx et al. 
2003; Mobasher et al. 2004; Wiklind et al. 2008). 
Franx et al. calculate a median age of 0.7-1 Gyr 
from a sample of J s — K s selected galaxies at pho- 
tometric redshifts ~ 3. A later paper by Wiklind 
et al. (2008) identifies galaxies at 3 > 5 using the 
3.6 fi m and K bands to detect the Balmer-break. 
They select 11 candidates with ages 0.2-1 Gyr, 
which implies extremely efficient star-formation 
processes in the early universe. Studies such as 
these provide evidence of evolved galaxies at high 
redshift, which is an important aspect to fully un- 
derstanding galaxy' evolution. 

This paper includes the use of optical to MIR 
observations based on sources extracted from the 
Hubble Ultra Deep Field Follow-up Near Infrared 
Camera and Multi-Object Spectrometer (HUDF- 
NICMOS) parallel observations (UDF05; Stiavelli 
2005), which is an extension of the NICMOS 
UDF Treasury Program (Thompson et al. 2005). 
Thompson et al. presented parallel observations 
of the HUDF with one subsection of the HUDF 
covered. Oesch et al. (2007, 2009) present results 
from this survey, showing observations of 3 > 5 
galaxies. These results are based on ACS obser- 
vations on the Thompson et al. (2005) NICMOS 
parallel fields. The results in this paper are galax- 
ies with redshifts 2 ~ 1.5 - 4, using the NICMOS 
observations, which are parallel to the ACS obser- 
vations in the Stiavelli (2005) UDF05 study. 

We describe the data reduction, source extrac- 
tion, error analysis and catalog selection in §2. In 
§3, the SED fitting method is explained in de- 
tail. In §4, We provide results for a selection 
of the data, which includes Monte Carlo sim- 
ulations to test the robustness of the SED re- 
sults. We have adopted the cosmological constants 
H 0 = 70 km s -1 Mpc -1 , f2 m = 0.3. and H A = 0.7 
throughout this research. 

2. The Sample 

The sample includes the multi- wavelength anal- 
ysis of galaxies observed in the NICMOS par- 
allel pointings in the Hubble Ultra Deep Field 
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Follow-up (UDF05; Stiavelli 2005). The UDF05 
project is an HST Large Program of ultra-deep 
ACS/WFC and NICMOS observations. As shown 
in Figure 1, the ACS/WFC observed the original 
HUDF parallel fields NICP12 and NICP34, im- 
aged by the NICMOS F110W (J 110 ) and F160W 
(Hieo) filters and are not used in this paper. The 
UDF05 also conducted parallel NICMOS observa- 
tions, where one of the pointings is on the orig- 
inal HUDF field (numbered regions in Figure 1). 
These parallel observations are the main source 
of photometry for our dataset. In addition to the 
multi-wavelength coverage, the reason for selecting 
the NICMOS parallel fields was to have coverage 
of multiple areas to reduce cosmic variance effects 
(Somerville et al. 2004). 

This section discusses the data reduction, 
source extraction, catalog construction, and data 
selection techniques that we used for these data. 

2.1. Data Reduction 

We extract sources from five parallel fields, in- 
cluding the HUDF, for the main catalog. Fig- 
ure 1 displays the fields for the HUDF and the 
HUDF parallels NICP12 and NICP34, along with 
the parallel NICMOS fields for the UDF05. Red 
boxes mark the 5 parallel UDF05 NICMOS fields. 
Each one has a reference number; for example, the 
HUDF pointing is Region 1. The total area cov- 
ered by the 5 regions is ~ 5.7 arcmin 2 . 

The NICMOS images in the F110W (Jno) and 
F160W (Hnso) filters were downloaded from the 
Multimision Archive at STScI (MAST) archival 
search tool. We followed the technique in Thomp- 
son et al. (2005) for the reduction of the images. 
The NIC3 camera has a pixel size of 0.2"pix _1 , 
which makes it impossible to resolve some of the 
more compact galaxies at higher redshifts. To in- 
crease the sampling for each pixel, we drizzled 
the multiple pointings in each field together us- 
ing the IRAF 1 tool MultiDrizzle to 0.09"pix _I . 
In the drizzling process, the MultiDrizzle param- 
eter PIXFRAC was set to 0.6 and SCALE to 
0.09/0.202863, which produced the 0.09" pix -1 . 
The number of combined images for each field, the 
calculated exposure times, and the drizzled image 
sizes are listed in Table 1 . The final processed im- 


1 Image Reduction and Analysis Facility is available at 

http://iraf.noao.edu/ 


ages are displayed in Figure 2. They are ordered 
by region number from top to bottom, and the left 
and right columns are the Jno- and Hmo-bands. 

This project consists of multi-wavelength anal- 
ysis, which includes observations in the optical 
(Advanced Camera for Surveys with the Wide 
Field Camera, ACS/WFC) and MIR (Spitzer 
Space Telescope with the Infrared Array Cam- 
era, SST/IRAC). For the optical wavelengths, 
we downloaded all data from MAST and used 
it as processed in the pipeline. For Region 1, 
the HUDF (Beckwith et al. 2006) ACS F606W 
( V 606 ) 7 F775W ( 2775 ), and F850LP ( 2 S5 o) images 
were acquired from the MAST archive. 

The other fields, Regions 2-5, are outside of 
the UDF and Great Observatory Origins Deep 
Survey (GOODS; Giavalisco et al. 2004) foot- 
print, so only the ACS F606W (Woe) and F850LP 
(zg 5 o) from the Galaxy Evolution From Morphol- 
ogy (GEMS; Rix et al. 2004) survey were avail- 
able. The GEMS footprint only partially covered 
Region 5. GEMS has a 5cr depth of mAB(Veo6) 
= 28.3 and mAB^sso) = 27.1. 

The ACS images (0.05"pix“ 1 ) were rebinned 
to the drizzled NICMOS pixel size (0.09"pix _1 ). 
They were then convolved with a NICMOS point- 
spread function (PSF) generated with the Tiny 
Tim (Krist 1993) modelling software. While this 
increases the error for the ACS photometry, it al- 
lows for better SED fitting since the magnitudes 
are measured over similar isophotal areas. 

The Spitzer/IRAC filters include four channels 
(3.6, 4.5, 5.8, and 8.0 pm). All data were down- 
loaded from the SST archive retrieval tool Leop- 
ard Version 9 2 . We obtained the Region 1 data 
from the GOODS Spitzer Legacy Science Program 
(Dickinson et al. 2004). The GOODS program 
reach 5 a depth at toab = 26.1, 25.7, 23.8 and 
23.7 for the 3.6, 4.5, 5.8, and 8.0 pm filters, re- 
spectively. We retrieved the images for Regions 
2-5 from the Extended CDF-S program by van 
Dokkum et al. (2005) (SIMPLE program). For 
all IRAC data, we analyzed the Post-BCD mosaic 
images that have pixel scales of 0.6" with no fur- 
ther processing. 


* Leopard is available at http://ssc.spitzer.caltech.edu/propkit/spot/ 
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2.2. Source Extraction 

We used SExtractor (SE Bertin & Arnouts 
1996) in dual-image mode for object detection. 
The primary image for NICMOS was the Hi6o fil- 
ter. Table 2 lists the parameters used in the SE 
input file. For ACS source extraction, the z 8 so 
image was the primary image. The GOODS SE 
detection parameters were used with a few adjust- 
ments for background. For IRAC source extrac- 
tion, the 3.6 //in image was used as the primary 
image and SE parameters came from Dickinson et 
al. (2004). 

For the NICMOS drizzled images, the edges 
have a much smaller integration time compared to 
the inner portions of the image. In order to mini- 
mize this, we created masks for the outer edges in 
each image. The drizzled weight maps do account 
for some of this effect (the exposure time largely 
determines the weight image values), but further 
masking was needed to reduce the problems asso- 
ciated with edge-detections. 

The photometric measurements are isophotal 
magnitudes for all ACS and NICMOS magnitudes, 
and 3" diameter apertures are used for IRAC 
sources. The IRAC detections are known for prob- 
lems with confusion due to its low resolution (i.e., 
overlap between sources that are too close to be 
resolved), so special care was introduced to han- 
dle this. We used a cross-catalog matching al- 
gorithm to indentify NICMOS objects that were 
multiple matches to IRAC detections and checked 
each multiple-match manually for blended sources. 
If a blended source was present, we masked (from 
SE segmentation maps), where possible, the ob- 
structing object and re-measured the magnitude. 
If it was not possible to separate the sources, we 
removed them from the catalog. This affected only 
a handful of objects. 

2.2.1. Photometric Errors 

For the NICMOS objects, the errors include 
those calculated from Poisson noise in the weight 
maps and systematic errors. The systematic errors 
were determined by simulating galaxies at random 
positions and comparing the SE measured magni- 
tudes with the input magnitudes (see §2.4). The 
error is then calculated by Am = (m, n — m out< ), 
where m m is the input (model) magnitude and 
m ou t is the SE measured magnitude. This is the 


same method applied by Giavalisco et al. (2004). 
The standard deviation of the Am distribution is 
~ 0.05 mAB, which is added in quadrature to the 
SE measured errors. We also added (in quadra- 
ture) 0.07 mAB for NICMOS as a lower limit for 
the error, based on the ACS systematic errors de- 
termined by Giavalisco et al. (2004). The reason 
is that NICMOS, which is much lower signal-to- 
noise and has less resolution, will introduce more 
error into the photometry. Therefore, the NIC- 
MOS measurements will have, as a lower limit, 
the ACS systematic error (T. Dahlen 2009, pri- 
vate communication). 

2.3. Catalog Construction 

To create the final catalog, all SE catalogs were 
matched by RA and DEC. The catalog is based on 
the NICMOS Hi6o-band detections. The multi- 
wavelength catalogs created from SE were cross- 
matched using an automated program, which in- 
cluded flags for multiple matches and cases with 
no matches. The multiple matched objects were 
checked manually and the catalog was adjusted 
accordingly. For example, in the case of confusion 
the competing objects were masked and the ob- 
ject of interest measured again. Also, objects de- 
blended in the ACS, may not have been deblended 
in NICMOS because the point spread functions are 
dissimilar in different bandpasses. In these cases, 
we re-measured the NICMOS magnitudes, so that 
SE would not deblend these. 

Any object missing in the Jno-band was re- 
moved, which proved to be a useful check, since 
many of the non-detections in the Jno-band were 
very faint and beyond the detection limits in the 
Hi6o-band image (mAB > 26 see §2.4). This re- 
sulted in a total of 630 objects. 

The NICMOS photometry was checked against 
Thompson et al. (2005) and we found an offset 
of ~ 2.3 for the Hi6o-band and ~ 1.5 for the Jno- 
band, where our data are brighter by this amount. 
This difference is probably due to: 1) background 
subtraction in the drizzling process; 2) different 
methods of estimating the background in source 
extraction. We added the offsets to calibrate the 
data with the Thompson et al. magnitudes. A 
check was also conducted by fitting SEDs to a few 
objects, which determined that the objects with 
the calibrated magnitudes resulted in more robust 
fitting. 
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2.4. Detection Completeness 

To calculate the completeness, we used the 
NIC3 Hieo-band filter images and inserted artifi- 
cially created galaxies into each region. We ex- 
tracted sources with SE based on the method 
above. We randomly generated 500 galaxies of de 
Vaucouleurs (disk galaxies with surface brightness 
distributions oc (r/r e ) 1/4 ) and elliptical profiles at 
varying luminosities and at random positions, us- 
ing the IRAF tool artdata. These were added to 
the Hj6o-band images. A polynomial was fit to 
the detection ratio and magnitude values. We re- 
cover 80% of the artificially generated sample at 
wtab = 25.4 and 50% at toab = 26.7. Figure 3 
is a plot of the detection ratio with tuab for the 
H 160 -band images. We calculated the complete- 
ness values at the 50% and 80% detection ratios 
by solving a polynomial fit (dashed line) at those 
points. 

Figure 4 displays the toab distribution for the 
V 606 , Z 850 > Jno, and Hmo filters of the actual de- 
tections (not artificially generated objects). The 
mAB are noticeably higher for the optical, espe- 
cially the Veoe-band, which is to be expected. One 
reason for this is that the signal-to-noise (S/N) is 
much higher for observations in the optical com- 
pared to the NIR. Also, depending on the rest- 
frame wavelength, a dusty or older population 
galaxy will register a higher flux for wavelengths 
red-ward of rest-frame wavelengths. For example, 
Figure 5 (discussed in §3) shows much brighter 
NIR rest-frame flux for the older 290 Myr and 1.4 
Gyr galaxies. For higher redshift galaxies (z > 1), 
this results in a larger magnitude in the V 606 -band 
compared to the Jno- and Hi 6 o-bands for the same 
object. 

3. Spectral Energy Distributions 

The SED modelling in this paper is based on 
Bruzual k Chariot (2003) (hereafter BC03) with 
Chariot k Bruzual (2007) (hereafter CB07) adap- 
tations for handling the TP-AGB phase (ther- 
mally pulsating asymptotic giant branch stars). 

In our data, the location of the Balmer-break 
(A RF = 3646 A) with respect to the observed 
wavelength is important, since the SED fitting 
method used for our photometry exploits this fea- 
ture to identify the best-fit SED. The Balmer- 
break is most prominent in A stars, since in O 


and B stars the hydrogen is mostly ionized, and is 
a useful age diagnostic for galaxies between 0.1 to 
1.0 Gyr. 

3.1. SED Parameters 

The parameters obtained from SED fitting rely 
upon proper redshift fitting. The method applied 
in this paper is a standard x 2 minimization tech- 
nique. Generally, the SED fitting is employed with 
a set of a priori conditions on the parameters de- 
termined by the spectra (metallicity, dust extinc- 
tion, starburst age, etc.). For this dataset, the 
observed photometric points from the catalog are 
compared to a set of template spectra generated 
from the BC03/CB07 libraries and with a large 
parameter space for redshift, starburst age (t s b), 
extinction (Eb-v), metallicity (Z) and star forma- 
tion e-folding decay time (r). While the parameter 
space allows for an unbiased and prior-free SED 
fitting, there are cutoffs which apply hidden priors. 
Also, the number of filters used for photometry is 
very important in determining the x 2 value, since 
the number of photometric data points limits the 
number of parameters that define the SED. 

The SED colors are evaluated in fixed observed 
bands (HST ACS: V 606 , * 775 , and z 850 ; HST NIC- 
MOS: Jno and Hi6o; Spitzer IRAC: 3.6, 4.5, 5.8, 
and 8.0 /urn). The ranges and cutoff's for the pa- 
rameters are similar to Wiklind et al. (2008): 
z = 0.1 — 9 with Az = 0.1, Eb-v = 0.0 — 0.95 
with AEb-v — 0.025, r = 0.0 — 1.0 Gyr with 
At = 0.1, and Z = 0.2, 0.4, 1.0, and 2.5 Z©. The 
starburst age ranges from t s b = 0.005 - 2.4 Gyr 
with steps of 5 Myr for f S b <0.1 Gyr and age steps 
of 100 Myr for t s b > 0.1. For comparison with the 
age of the Universe, at z ~ 2.7 the Universe was 
approximately 2.4 Gyr. 

To account for inter-galactic medium (IGM) 
attenuation due to neutral hydrogen absorption, 
the Madau (1995) prescription is used for the 
mean IGM opacity. Madau provides transmis- 
sion curves for high-; (2.5 < z < 4.0) galaxies as 
a function of observed wavelength, which can be 
used to correct the flux based on the redshift or 
rest-frame wavelength. The magnitudes can in- 
crease by ~ 2 in the redshift range 2.5 < z < 4.0. 
The dust extinction Eb-v was calculated using 
the attenuation law from Calzetti et al. (2000) 
and is parameterized through the Eb-v color in- 
dex. The bolometric luminosity (Lboi) is derived 
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by integrating over the entire wavelength range 
(91 A to 160 i-im). The SED output provides Lb 0 i 
with and without dust corrections, which gives a 
method for estimating the FIR luminosity (Lffr). 
Calzetti et al. (2000) define the color excess based 
on comparisons between the UV and FIR flux for 
starburst galaxies. They determine that ~ 30% 
of the bolometric flux is from the UV-NIR wave- 
length range and the rest is emitted in the FIR. By 
assuming that dust absorption of the UV-optical 
is re-emitted in the FIR, the FIR luminosity can 
be approximated by calculating the excess. All 
SED values for FIR luminosity in this paper are 
given by L F ir = L Bo i - L Bo i e , where L Bo i e is the 
bolometric luminosity corrected for dust. 

The SSPs in BC03 give the stellar mass-to-light 
ratio M*/Ly based on the V-K and B-V colors 
and metallicity. As reddening and metallicity in- 
creases, so does the M*/Ly. The physical rea- 
son for this is that, for a fixed stellar mass, lower- 
ing the metallicity causes stars to evolve at higher 
effective temperatures and luminosities (Schaller 
et al. 1992; Fagotto et al. 1994; Girardi et 
al. 2000). After identifying the best-fit model, 
the stellar mass M* is calculated by simply mul- 
tiplying L B oi by M*/Ly. The star formation rate 
(SFR) is calculated through an exponential declin- 
ing function, based on the e-folding time (r) and 
the starburst age (t s b) and results are discussed in 
§4.2. 

Several million models are pre-calculated, span- 
ning the predefined parameter spaces. The result- 
ing SEDs are integrated through filter response 
functions for the appropriate filter, which are then 
best-fit with the observed photometry using y 2 
minimization. 

4. Results 

We created a subsample of 301 objects out of 
630 in the main catalog. This selection is based 
on the number of free parameters used for fitting 
the SEDs (described in §3). Since there are five 
parameters, at least five photometric values must 
be present for a proper fit, which resulted in 301 
objects. Most of these objects had no IRAC detec- 
tions, which makes the sample limited to the IRAC 
sensitivities. Because the sources axe selected from 
the Hi6o-band, the sample is more sensitive to stel- 
lar mass than UV luminosities. 


4.1. Validity of Photometric Redshifts 

Figure 6 displays the distribution of photomet- 
ric redshifts, or z p h 0 t , obtained from the SED fit- 
ting. There are noticeable peaks that indicate 
both the strengths and weaknesses of the Balmer- 
break technique. Because the Balmer-break oc- 
curs at 3646 A, for any galaxy at z < 0.5 (Xrf < 
6000 A) the break will not be detected in this sam- 
ple. This is a selection effect from the limited 
optical photometric bands; the ACS-V 606 cov- 
ers the edge of this break at this redshift limit. 
Without the ACS-B 435 -band the break cannot be 
determined at lower redshifts (see Figure 5). A 
strength for the technique with this sample occurs 
for 2 > 1.3 (A rest > 9200 A), where the ACS - 2350 
and NICMOS coverage is applicable for 1.3 < 2 < 
3.5. At even higher redshifts, the IRAC images 
provide the signatures of the Balmer-break, but 
significant problems may arise from the large er- 
rors and confusion intrinsic with infrared imaging, 
due to the low resolution of Spitzer. 

In order to determine the efficacy of the 
2 p hot results, we matched multiple spectroscopic 
redshift catalogs with the main SED catalog and 
found 16 objects (Objects 7, 13, 15, 21, 24, 25, 38, 
58, 71, 76, 78, 102, 119, 133, 144, and 204). Fig- 
ure 7 shows the comparison between the z p hot and 
spectroscopic redshifts (z spec ) for these 16 ob- 
jects. They are discussed in more detail in the 
following section. Figure 7 compares z spt . c and 
^phot using A — (2 p hot ^spec)/(l ± z spec ). This 
is a standard method to determine the fractional 
error (see Mobasher et al. 2004). The figure lists 
the median error, < A >, the rms scatter a. The 
small number of objects biases our results, but 
there is a clear separation between z p hot obtained 
for 2 < 1 and 2 > 1. The objects at higher red- 
shifts (2 > 1) have z P hot that agree well with the 
spectroscopically confirmed redshifts. Out of the 
16 galaxies, 75% have A within ±0.2. 

4.2. Objects with Spectroscopic Redshifts 

Figures 8-9 show eight greyscale images for the 
16 objects that have available spectroscopic red- 
shifts. The images in the top row are: ACS Veo6, 
zgso, and NICMOS Jno, Hi6o- The IRAC images 
are in the bottom row for Chi , Ch2, Ch3, and 
Ch4 (3.6, 4.5, 5.8, and 8.0 pm) as labeled. A 1" 
ruler marker is provided for each image. 
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The best-fit SEDs are shown in the Figure 10. 
The SED shows the observed and rest-frame wave- 
lengths marked as the bottom and top abscissa, 
respectively. The magnitudes are listed in AB, 
which is the same as the photometric data in the 
main catalog. The two different SED lines refer to 
the flux corrected for dust and flux without dust 
corrections (red and blue, respectively in online 
color figure). Plots with only one line have Eb-v 
= 0, so no dust corrections are calculated. This is 
the case for objects 38, 58, 76, 144, and 204. 

Table 3 gives a list of the 16 objects with 
spectroscopic redshifts and their best-fit param- 
eters. Of special note are the 11/16 galaxies with 
A < 0.2, which are specified with an asterisk. 
Galaxies with the highest redshifts are Objects 24 
and 38 at z spec = 3.6 and 5.95. 

In Figure 8, Object 24 has a very bright lentic- 
ular structure with a young (f s t, = lOMyr) dusty- 
starburst SED template. It appears as an edge-on 
disk in the ACS V 606 and z 850 , and has a promi- 
nent bulge in the NICMOS images. The ACS 
corresponds with the UV rest-frame, while the 
NICMOS corresponds with the optical rest-frame. 
This object has an effective radius of 0.54" in the 
observed NIR, or ~ 4.5 kpc at z = 3.6. It is ex- 
tremely massive, with a stellar mass of ~ 10 12 Mg. 
Based on the SED-derived Lfir, this is a ULIRG 
at 2 ~ 4. A Spitzer MIPS detection is incon- 
clusive, because the pixel scale is 2.45"pix -1 ; one 
pixel is twice the size of the object. It is the most 
massive and FIR-luminous out of the 11 spectro- 
scopic redshift objects. 

Object 38 is not well resolved, but appears to 
be a clumpy, compact galaxy. In Figure 8 the de- 
tection is very faint in the optical-NIR and “drops 
out” in the Vgoe-band. It also has SED features 
of a young (t s b = lOMyr) starburst, but with no 
dust correction. The effective radius is 0.27" in 
the observed NIR, which is equal to ~ 1.5 kpc at 
z — 5.9. It has a stellar mass of ~ 10 s Mg and 
Lfir ~ 10 10 Lq. This is a compact dust-free star- 
burst galaxy and fits all the properties of being an 
LBG at z ~ 6. 

Objects 76, 102, 119 and 133 have z spec ~ 1.3. 
All but one (Object 76) have IRAC photometry. 
In Figure 9, Object 76 shows itself as a very thin, 
wispy feature in the Vgog and zgso (rest-frame 
FUV and U-band). In the NICMOS images, the 
galaxy pops out to almost a magnitude brighter 


(rest-frame optical). The effective radius is 0.62" 
in the observed NIR, or ~ 5.23 kpc at z = 1.3. 
The SED displays a dust-free older star-forming 
galaxy (t S b = 100 Myr), where the Balmer-break is 
stronger from the O and B stars moving out of the 
main sequence. It has a stellar mass of ~ 10® M@ 
and Lfir ~ 10 9 Lg. 

Object 102 (Figure 9) also has a thin appear- 
ance in the ACS, or rest-frame UV, images. In 
the NICMOS, or rest-frame optical, an outer disk 
of older population stars appears. The effec- 
tive radius is 0.42" in the observed NIR, or ~ 
3.5 kpc at z = 1.3. The SED is modeled by a 
very dusty ( Eb-v — 0.375), highly star-forming 
(SFR a „ e ~ 49M 0 yr -1 ) plus an older starburst age 
(t S b = lOOMyr). The Balmer-break is a tiny decre- 
ment in the zgso-J i 10 observations. It has a stellar 
mass of ~ 10 9 Mq and Lfir ~ 10 n Lg. 

In Figure 9, Object 119 appears bright in the 
rest-frame UV as an irregularly shaped, inner core 
in the ACS filters. Its reddened outer shell appears 
bright in the NICMOS images. There is confu- 
sion in IRAC with a bright feature in the lower 
left. This galaxy is compact with an effective ra- 
dius of 0.23" in the observed NIR, or ~ 1.9 kpc at 
z = 1.3. The SED is fit to a young dusty star- 
burst (t sb = 30 Myr) and displays a bright MIR 
bump, which may be falsely augmented by confu- 
sion. The dust extinction ( Eb-v — 0.375), Lfir 
(~ 10 u L 0 ) and M* (~ 10 lo Mg) are all similar 
to Object 102, which also has a slight MIR bump. 

For objects 15, 25, 58, 71, and 78 the SED fits 
are remarkably close and redshifts are at z ~ 1 
(see Figure 10). Object 15 is clumpy and dis- 
turbed morphologically. It appears to have been 
in a recent merger in the ACS images, where very 
patchy clumps of star-forming regions are seen in 
the rest-frame NUV (observed V 6 oe)- The galaxy 
becomes a bulge-dominated elliptical in the rest- 
frame optical (observed Higo)- Object 71 also 
has a clumpy rest-frame NUV morphology and 
becomes an extremely bright rest-frame optical 
through MIR source (observe zgso-IRAC 4). The 
Balmer-break has blended with the 4000 A break, 
which fits an older 1.02 Gyr galaxy with little dust 
extinction. The remaining objects have very sim- 
ilar morphologies at different wavelengths. 

All objects with insufficient fits for z P hot values 
occur for z spec < 1.3, which as discussed above 
in §3 is a weak region for the Balmer-break to be 
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detected with this photometric dataset. Addition- 
ally, many of the IRAC detections have confusion 
or are non-existent, which causes a crude fit. Ob- 
ject 7 is an example, where no IRAC detection is 
present and z p hot = 4.3 for a z s pcc = 103 object. It 
should also be noted that objects 7-133 with spec- 
troscopic redshifts are in Region 1, which is part 
of the HUDF. We obtained i 775 -band photometry 
for this pointing, which extended the SED cover- 
age between the Veoe- and £ 850 -bands. 

Three objects (21, 24 and 78) have unrealistic 
star formation rates (see Table 3). They are fit 
with SEDs that have dusty, instantaneous bursts 
(r = 0) with starburst ages 10 Myr. This is 
tested in the Monte Carlo simulations in the fol- 
lowing section. By fixing z pho t > thereby reducing 
the number of free parameters, a more likely age 
and star formation rate is generated. Even though 
Zphot is correct for object 24, the age and star for- 
mation rate is unlikely to be real. The results for 
Zphot and mass are more robust than the age and 
star formation rate, due to SED fitting methods. 
Lee et al. (2009) determine that the stellar mass 
is underestimated by 19%-25%, because of older 
stellar populations being hidden by recent star for- 
mation. Furthemore, they find the star formation 
rate is less robust due to various star formation 
histories in galaxy templates, and the degeneracy 
between age and dust extinction. 

4.3. Monte Carlo Simulations 

The foremost potential source of systematic er- 
rors in SED fitting is the reliability of the pho- 
tometric redshifts, which we tested, as mentioned 
above, by comparing with the available spectro- 
scopic redshifts. Indirect testing may also be per- 
formed for all the objects through Monte Carlo 
simulations. These simulations also check for the 
reliability of other parameters that are also vulner- 
able to erroneous values. The Monte Carlo simu- 
lations provide the probability distribution of the 
parameters. The simulations consist of pseudo- 
random re-sampled magnitudes, within the given 
photometric errors. Slight shifts in the magnitudes 
can create different stellar population models to 
be fit to the input SEDs, but should give a normal 
distribution around a peak for a robust fit. 

Table 4 shows the parameter values for the 
Monte Carlo simulations for all objects with spec- 
troscopically confirmed redshifts. For each ob- 


ject, there are two sets of parameters. The first 
set allows for the redshift to be a free param- 
eter in the fitting (top). The second set gives 
results for a constant redshift at z spcc (bottom). 
These results include 500 simulations for each ob- 
ject, giving the probability density of a certain 
value. All values presented are the distribution 
peak value unless indicated otherwise. Columns 2 
and 3 list the redshifts for the Monte Carlo distri- 
bution peak (z mc ) and the best-fit (z p hot)- Column 
4 gives the percentage of the simulated fits that 
obtained the z p hot value. Columns 5-10 are: the 
Eb-v , the metallicity, Lb 0 i without dust correc- 
tion, M*, mean SFR ave , mean t s t,. Note that the 
Zphot obtained from the SED fitting does not nec- 
essarily agree with the most likely median value 
in the Monte Carlo simulations. This is probably 
caused by the relatively flat SED and the number 
of degrees of freedom in the SED fitting process 
compared to the number of photometric points. 
Setting the redshift to a constant shows more re- 
alistic values in terms of the star formation rates, 
indicating that the number of photometric points 
is a cause of the multiple peaks in the Monte Carlo 
simulations. 

4-3.1. Three Study Cases 

The parameter distributions in Figure 11 often 
show multiple peaks. The presence of the multi- 
ple peaks indicates that there is a certain proba- 
bility for the SED to fit that value. We discuss 
the results for three galaxies based on the several 
possible parameters that could 

Object 24 (top panel in Figure 11) is an ex- 
ample of an object where the Monte Carlo distri- 
bution peak redshift, z mc , agrees well with both 
Zphot and z spec • The redshift is narrowly dis- 
tributed, while the double peaks for other param- 
eters indicate a significant likelihood for the two 
values. The stellar mass, dust extinction, and lu- 
minosity result in best-fit values within the dis- 
tribution of these peaks and are not too different 
from the maximum from the Monte Carlo distri- 
bution. 

Object 25 (middle panel in Figure 11) is one of 
the typical results, where z mc (3.8) does not agree 
well with Zphot (1-1) and z spec (1.09). All the pa- 
rameters have multiple peaks, where the best-fit 
result is located within one of the peak distribu- 
tions (not necessarily the most populated). 
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Object 144 (bottom panel in Figure 11) is an 
example of one of the worst fits, where z mc (5.7), 
Zphot (0.1) and z spC c (0.54) do not agree and the 
parameters have very noisy distributions. The 8.0 
/mi IRAC magnitude has a very large error and 
this is very likely what contributes to the compli- 
cated distribution for redshift, which then makes 
the remaining parameters difficult to fit. 

As a representative sample, the spectroscopi- 
cally confirmed objects have 4 good fits similar 
to Object 24, 8 typical fits similar to Object 25 
and 4 bad fits similar to Object 144. Statis- 
tically this means that the SED fitting method 
works for approximately 75% of the objects at 
z > 1.3. A larger number of objects with spectro- 
scopically confirmed redshifts is needed to validate 
this trend. However, the Monte Carlo method is a 
powerful tool in evaluating the SED fits and pro- 
vides confidence when evaluating the results of the 
fit. 

4.4. Galaxies at z~1.3 and 4 

We selected galaxies, from the SED output cat- 
alog, that are in the redshift ranges 1.3 < z <2 
(35 objects) and 3 < z < 4 (31 objects) and that 
also had at least one IRAC detection. The pur- 
pose of the redshift selection is that these red- 
shifts have been shown (see §4.1) to be appropriate 
for Balmer-break detection with our photometric 
data. The IRAC photometry helps constrain the 
fit for the rest-frame NIR to MIR stellar popula- 
tion and tend to have better agreement with z spe c - 
Six of these objects (Objects 38, 58, 76, 102, 119, 
and 133) have z spec and agree with the trends de- 
scribed below. 

Figure 12 shows six histograms with distribu- 
tions of SED parameters. The diagonally hatched 
bars depict galaxies with spectroscopic redshifts at 
1.3 < « < 2 (objects 76, 102, 119, 133) and z > 3 
(objects 24 and 38), respectively. Several notice- 
able trends are immediately apparent. In the top 
right plot, it is clear that the higher redshift galax- 
ies have less dust extinction. This is a selection 
effect, because the rest-frame UV is observed by 
optical wavelengths for the higher redshift popu- 
lation. The UV light is highly absorbed by dust. 
Any object too faint from dust extinction, and not 
detected in the UV rest -frame would not have the 
required photometric points for SED fitting. One 
way to solve this is by having more IR photometry 


and JWST will provide the photometry necessary 
to extend the SED at these redshifts. 

The luminosities and star formation rates show 
brighter sources for higher redshifts. The increase 
in luminosity at z ~ 3 — 4 agrees with Bouwens 
et al. (2007), which show a peak SFR density 
at this redshift range. Although, one might argue 
that this is Malmquist bias, it has been shown that 
the luminosity increase for z ~ 2 — 6 galaxies is 
not a selection effect (Bouwens et al. 2004, 2007). 
These authors used the sizes and magnitudes of 
galaxies in the UDF and GOODS and found that 
the GOODS objects represented the same sample, 
even though it is not as deep of a survey as the 
UDF. 

The stellar mass appears to be hovering around 
the same median point for both redshift ranges. 
In general, the majority of bolometric luminos- 
ity originates from massive stars. From where is 
the higher luminosity originating at higher red- 
shifts? This is most likely caused by the limits in 
the Salpeter IMF used in the SED code. In gen- 
eral, the majority of light measured from galaxies 
is produced from a few massive stars, while the less 
massive stars contribute most of the stellar mass. 
Furthermore, previous analyses have shown that 
there is a shortfall in ionizing photons to produce 
Ly-a flux at the re-ionization epoch (z ~ 7) using 
current stellar population models (e.g., Bromm, 
Kudritzki, k. Loeb 2001; Schaerer 2003; Bunker 
et al. 2009). One solution proposed includes in- 
creasing the IMF mass limit to a top-heavy IMF 
or to decrease metallicities. However, Franx et al. 
(2008) show that the stellar mass density increases 
with redshift for a given mass. Also, it may be 
due to evolution of the fundamental plane, where 
the mass-to-light ratio decreases with redshift and 
can be modelled by the evolution of a simple stel- 
lar population (see van Dokkum & van der Marel 
2007, and references therein). 

5. Summary 

We extracted sources observed in the paral- 
lel NICMOS UDF05 fields and selected galaxies 
to calculate SEDs. A total of 630 galaxies were 
cross-matched using the ACS filters Vgoe arid zsso, 
NICMOS filters Jao and Hmo, and Spitzer IRAC 
channels 1-4. The photometry are sensitive to 
stellar mass, because the sources are H^o-band 
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selected A selection of 301 objects were modeled 
with spectral energy distributions using x 2 min- 
imization. Photometric redshifts, dust extinc- 
tion, stellar masses, bolometric luminosity, star- 
burst age and metallicity were estimated to ana- 
lyze galaxy evolution over a select range of red- 
shifts. 

The main results are: 

1. The redshift distribution for 301 photomet- 
ric redshifts shows a patchy, non-normal dis- 
tribution. This is a result of the wavelength 
coverage, where the data is not sensitive to 
strong spectral features, such as the Balmer- 
break, which is sensitive to stellar age. It 
was shown with simple stellar population 
models that the break might be detected at 
z > 1.3 with the given photometry. 

2. Comparisons with 16 spectroscopic redshifts 
give 75% within A < 0.2. These photomet- 
ric redshifts are in accord with the sensitiv- 
ities expected from the Balmer-break with 
this dataset. The accuracy of the photomet- 
ric redshifts may be improved with color se- 
lection criteria. 

3. A case study of the Monte Carlo simula- 
tions determines that the parameters from 
the best-fit SEDs are robust for the redshift 
ranges 2 > 1.3. 

4. We analyze a subsample of galaxies at pho- 
tometric redshifts 1.3-2 (35 objects) and 3- 
4 (31 objects) and determine that Lb 0 i and 
SFR increase significantly from z ~ 1.5 to 
4. This agrees with previous results by 
Bouwens et al. (2007), and is not a selec- 
tion effect. Dust extinction decreases, but 
this is due to a selection effect, where the 
rest-frame UV and optical are observed in 
the optical and NIR at z ~ 4. The dust will 
absorb these wavelengths, which limits the 
photometric detections. 

5. Objects 24 and 38 are the highest redshifts 
that are spectroscopically confirmed at z ~ 
3.7 and 6. The SED model for Object 24 
gives Lfir = 10 13 L©, which would make it 
a ULIRG at z = 3.7. 

6. Objects 76. 102, 119, and 133 are objects 
with spectroscopically confirmed redshifts at 


z ~ 1.5. Objects 102, 119, and 133 have 
starburst ages f s b < 50 Myr, and are ac- 
tively forming stars. Object 76 is probably 
an older 100 Myr starburst with little active 
star formation. 

In future surveys, Objects 133, 24 and 38 are 
good candidates for JWST and ALMA. Object 
133 is a galaxy that formed when the Universe was 
~ 4.5 Gyr old. An instantaneous starburst with 
very high metallicity (2.5 Z O dot) fits the SED, im- 
plying previous multiple starbursts. JWST would 
easily resolve the rest-frame MIR features which 
are not clear in IRAC, due to confusion and low 
resolution. Similarly, Object 24 is an example of a 
galaxy that formed when the Universe was ~ 1.6 
Gyr old. It is possibly a ULIRG and is an example 
of a progenitor to nearby massive ellipticals. Ob- 
ject 38 is a more evolved galaxy (~ 1 Gyr) at z ~ 
5.8 and is possibly among the first generation of 
galaxies. A follow-up study on these galaxies and 
similar types, with the JWST and ALMA observa- 
tories would provide better resolution in the NIR- 
millimeter wavelengths (optical-FIR rest-frames), 
which is important for distinguishing dusty star- 
forming galaxies, dust-obscured AGN and more 
evolved stellar populations. The new observato- 
ries will reveal the nature of these distant objects 
and how they connect with the local universe. 

It is important to include high redshift observa- 
tions of the more evolved stellar populations, be- 
cause there has been a heavier emphasis on detect- 
ing LBGs and dusty starbursts at high redshifts, 
where more modest star-forming galaxies are not 
selected. This biases galaxy evolution models to- 
ward a younger population of galaxies. A com- 
parison of the young, starburst galaxies with the 
more modestly star-forming galaxies can lead to 
an extensive view of galaxy evolution. The older 
stellar populations at high redshifts, after all, were 
the new stars of the first galaxies. 
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Fig. 1. — Footprints of the UDF05 NICMOS parallel fields. The image displays the original HUDF field 
with the parallel fields NICP12, NICP34, and NICP5 overlayed on a DSS image (Stiavelli 2005). The small 
boxes, numbered 1-5, indicate the UDF05 parallel fields observed while ACS was observing the NICP fields. 
The numbers are referenced throughout the paper as Regions 1-5. 
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Fig. 2. — UDF05 NICMOS detection images. The ten images display the five UDF05 NICMOS parallel fields 
after data reduction and drizzling processes. The left panels show the F110W filter and the right panels 
show the F160W filter. Regions 1-5 are labeled for each set. 


Completeness for F160W 



Fig. 3 — Detection Completeness for NICMOS F160W. The data are the detection ratio at magnitude 
(M A b) for the NICMOS Hi 60 -band along with a polynomial fit (dashed line). Completeness at 80% and 
50% are listed in the upper right corner and are also marked with the vertical lines. 



Fig. 4 .— Magnitude histograms. The M AB distributions for the real data sample in the ACS V 60 6 and z 850 
filters (top), and for the two NICMOS Jn 0 and Hi 00 filters (bottom). 
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Fig. 5. — Simple stellar population templates for a galaxy with solar metallicity and evolved to three different 
ages: 5 Myr (blue solid line), 290 Myr (green dash-dotted line), and 1.4 Gyr (red dashed line). The four 
plots show the spectral lines redshifted for 2 = 0.5, 1.3, 3.0,and4.0 as labeled. Rest-frame wavelengths (Xrf) 
are labeled along the top axis, corresponding with the labeled observed wavelength (Ao) along the bottom 
axis; note that the wavelengths are logscale to embellish the spectral features, particularly the Balmer-break 
feature, for clarity. The bandwidths for the ACS V 6 oe , zsso ■ NICMOS Jno, Hmo, and IRAC chi, ch2, ch3, 
and ch4 filters are shown in the shaded regions. Note the redshifting of the Balmer-break at = 3646 
A (0.3646pm). The templates are part of the BC03 model libraries and have exponentially decaying star 
formation. 



Fig. 6. — Distribution of photometric redshifts. The z p h 0 t based on the results of the SED fitting for all 
301 objects (black line), and the 231 objects with IRAC photometry (grey filled) show peaks at particular 
redshifts. 
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Fig. 7. — Comparison of z p hot and 2 spec . The data show the fractional error A as a function of 2 spec - The 
rms scatter a and average error < A >= (z p hot — z Bpec )/(l + z sp ec ) are listed in the top right comer. We 
determine that objects at 2 > 1.3 are much more likely to be correctly measured than at lower redshifts, 
because of the bandwidth coverage (or lack thereof) in the dataset (see Figure 5). 
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Fig. 8. — The 8 postage stamp panels for spectroscopically confirmed objects display the separate filter 
images for the V 60e , 2 850 , Jno, H 160 , and IRAC channels 1-4 as labeled. Note the 1" ruler on each image 
for scale comparisons. 
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Fig. 10. — The best-fit SEDs for spectroscopically confirmed objects. Included in each plot are the photo- 
metric and spectroscopic redshift values. See Table 3 for the calculated SED parameters. The blue/red line 
is the flux with/without dust (see online version for color plots). 


20 





LOG(M*/M.) E(B-V) LOG(L Ilol ./L . ) 


Obj 144 atsv 

♦ 

* 

N1CH 

4m 

LOG(<SFR>) (M yr 1 J 




LOG(M*/M . ) E(B-V) LOGSL lMi ,/L . ) 


Fig. 11. — Monte Carlo simulation results for Objects 24. 25, and 144. The five Histograms show different 
SED parameters determined from random resampling through Monte Carlo simulations. The parameters 
are from left to right: redshift, average star formation rate, log(M*), Eb-v, and log(Leoi) (no dust). The 
dotted line shows the results based a constant 5 p hot ^ 















Fig. 12. — SED parameters for redshift-selected galaxies. From left to right, Eb-v > Lb 0 i (without and with 
dust extinction), Lfir, average SFR and M* are plotted for galaxies in the redshift ranges: 1.3 < z < 2 
(35 objects) and 3 < 2 < 4 (31 objects). The diagonally hatched bars indicate galaxies with spectroscopic 
redshifts at 1.2 < 2 < 2 (objects 76, 102, 119, and 133) and 2 > 3 (objects 24 and 38). 
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Table 1 

Exposure Properties of the UDF05 NICMOS Parallel Fields 


Region 

^exp 

sec 

# Images 

Area 

arcmin 2 

1 

52158 

40 

0.95 

2 

17919 

14 

1.32 

3 

15359 

12 

1.23 

4 

23039 

18 

1.72 

5 

60157 

47 

0.86 


Table 2 

Source Extraction Parameters Used for the NICMOS catalog 


Parameter 

Value 

DETECT -MINAREA 

16 

THRESH-TYPE 

RELATIVE 

ANALYSIS.THRESH 

1.2 

DEBLEND-NTHRESH 

32 

DEBLEND JvIINCONT 

.005 

WEIGHT-TYPE 

WEIGHT-MAP 

WEIGHT-IMAGE 

mask. fits, weightmap.fits 

FILTER-NAME 

gauss-1. 5-3x3. conv 
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Table 3 

Best Fit SED Parameters for Objects with Spectroscopic Redshifts 


Object 

2»p«C 


E B -v 

tSB 

Gyr 

r 

Gyr 

z 

Z 0 

log(L Bo |) 

Lq 

log(L F m) 

L© 

log(M*) 

M e 

M/L 
M G L- 1 

SFRi “ 
M 0 yr -1 

SFR 2 a 
Mg yr -1 

SFR. ave 
M-, yr” 1 

X 1 

r* 

1.03 

4.3 

0.075 

0.10 

0.8 

0.2 

11.05 

11.25 

933 

0.007 

20 

22 

21 

4.05 

13- 

0.38 

1.1 

0.225 

0.05 

0.0 

2.5 

11.22 

11.47 

10.27 

0.040 

286 

472 

371 

2.40 

15*' 1 

1.09 

1.1 

0.3 

0.29 

0.0 

1 

10.83 

10.85 

10.39 

0.178 

15 

261 

86 

3.41 

21 c 

0.76 

4.0 

0.275 

0.01 

0.0 

0.2 

12.26 

13.31 

12.02 

0.048 

> 5000 

> 5000 

> 5000 

3.62 


3.67 

3.7 

0.3 

0.01 

0.0 

0.2 

11.88 

13.00 

11.71 

0.048 

> 5000 

> 5000 

> 5000 

5.21 

25* d 

1.09 

1.1 

0.05 

0.10 

0.0 

2.5 

10.28 

9.58 

9.23 

0.075 

10 

27 

17 

3.74 

38* d 

5.95 

5.8 

0.0 

0.01 

1.0 

0.2 

10.44 

10.58 

7.82 

0.001 

13 

13 

13 

18.48 

58*' 1 

0.77 

1.1 

0.0 

0.05 

0.0 

2.5 

9.89 

7.09 

8.49 

0.040 

5 

8 

6 

21.46 

71* b 

0.95 

1.0 

0.025 

1.02 

0.0 

1 

11.12 

9.86 

10.75 

0.408 

0 

565 

56 

5.51 

76* d 

1.30 

1.1 

0.0 

0.10 

0.0 

0.4 

9.42 

8.96 

8.43 

0.075 

2 

4 

3 

18.78 

78* b 

0.62 

0.9 

0.325 

0.10 

0.0 

0.2 

10.85 

11.38 

11.18 

0.489 

875 

2403 

1512 

7.75 

102* b 

1.27 

1.1 

0.375 

0.01 

1.0 

2.5 

10.39 

11.51 

8.68 

0.001 

48 

49 

48 

3.05 

119* d 

1.22 

1.2 

0.375 

0.03 

0.0 

2.5 

10.38 

10.94 

9.45 

0.025 

80 

108 

94 

12.43 

133" J 

1.31 

1.2 

0.025 

0.05 

0.0 

2.5 

10.14 

9.33 

8.81 

0.040 

10 

16 

13 

3.50 

141 b 

0.54 

0.1 

0.0 

0.72 

0.2 

0.2 

8.07 

6.68 

7.07 

0.098 

0.002 

0.06 

0.02 

1.70 

204 h 

0.18 

4.2 

0.0 

0.01 

1.0 

2.5 

12.27 

11.88 

9.57 

0.001 

366 

369 

368 

9.67 


'These objects have |A| < 0.2 

"SFRi is calculated at the age t = t„b! SFRj is calculated at the age t = 0. 
b z»pec were obtained from the WDS CDF-S (Le Fevre ft al. 2004) 
c 2»p*e were obtained from Szokoly et al. (2004) 

d - 2 «pe<- were obtained from the GOODS VLT FORS2 survey (Vanzella et al. 2006). 



Table 4 

Parameter Values from Monte Carlo Simulations 


Object 

r a 

‘me 

, b 

2 phot 

Percentage c 

2 phot 

< Eb-v > d 

z 

( Z ©) 

L ] 3 o 1 

(L 0 ) 

M* 

(M 0 ) 

< SFR ave > d 
( Mq yr" 1 ) 

< t« b > d 
(Gyr) 

7 

4.2 

4.3 

12.0 

0.04 

0.2 

11.2 

8.2 

103 

0.42 


1.0 

- 

- 

0.145 

0.2 

7.8 

7.6 

1 

1.93 

13 

4.0 

1.1 

0.0 

0.27 

1 

14.1 

11.1 

> 5000 

0.06 


0.3 

- 

- 

0.821 

0.2 

9.2 

8.9 

1 

2.04 

15 

1.0 

1.1 

11.2 

0.36 

2.5 

10.8 

10.3 

> 5000 

0.59 


1.0 

- 


0.901 

0.2 

9.6 

9.4 

6 

1.03 

21 

4.1 

4.0 

24.2 

0.19 

0.2 

13.4 

10.4 

> 5000 

0.04 


0.8 

- 

- 

0.551 

0.2 

9.5 

9.1 

1 

2.47 

24 

4.0 

3.7 

9.0 

0.16 

0.2 

11.8 

10.3 

> 5000 

0.05 


3.6 

- 

- 

0.461 

0.2 

9.1 

8.9 

1 

2.03 

25 

3.8 

1.1 

3.6 

0.27 

0.2 

13.0 

10.1 

> 5000 

0.15 


1.0 

- 

- 

0.518 

0.2 

8.9 

8.7 

1.2 

0.72 

38 

5.7 

5.8 

16.4 

0.01 

0.2 

10.7 

7.7 

18 

0.02 


5.8 

- 

- 

0.402 

0.2 

6.5 

6.4 

0.03 

0.52 

58 

4.0 

1.1 

1.6 

0.06 

0.2 

11.2 

9.3 

272 

0.13 


0.8 

~ 

~ 

0.144 

0.2 

8.0 

7.8 

0.10 

1.61 

71 

0.8 

1.0 

4.8 

0.21 

1 

10.9 

10.6 

> 5000 

1.2 


1.0 

- 

- 

0.818 

0.2 

9.9 

9.6 

3 

2.38 

76 

4.0 

1.1 

5.8 

0.07 

2.5 

10.9 

8.1 

65 

0.12 


1.2 

- 

- 

0.181 

0.2 

7.0 

6.9 

0.2 

1.09 

78 

4.4 

0.9 

1.0 

0.26 

1 

13.7 

10.7 

> 5000 

0.46 


O.G 

- 

- 

0.609 

0.4 

9.6 

9.3 

1 

2.87 

102 

3.5 

1.1 

0.6 

0.18 

0.2 

11.7 

11.3 

> 5000 

0.20 


1.2 

- 

- 

0.639 

0.2 

8.4 

8.2 

0.2 

0.82 

119 

3.6 

1.2 

6.6 

0.24 

0.2 

11.7 

10.8 

> 5000 

0.46 


1.2 

- 

- 

0.693 

0.2 

8.4 

8.2 

0.3 

0.52 

133 

0.1 

1.2 

3.0 

0.07 

1 

11.2 

9.4 

498 

0.09 


1.2 

- 

- 

0.207 

0.2 

7.6 

7.3 

0.1 

0.85 

144 

5.7 

0.1 

10.5 

0.04 

0.4 

12.5 

9.3 

532 

0.02 


0.5 

- 

- 

0.449 

2.5 

8.2 

6.2 

0.2 

0.77 

204 

6.0 

4.2 

0.0 

0.02 

2.5 

12.4 

10.4 

1805 

0.02 


0.2 

- 

- 

0.833 

2.5 

9.3 

7.1 

1 

0.69 


a 2 mc is the peak value of the redshift in the Monte Carlo simulations. 
b ~phot is the best-fit redshift measured from SED fitting. 
c The percentage of Monte Carlo simulations that resulted in the 2 P hot value. 
d < Eb-v >, < SFRave > and < t B b > are mean values. 
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